Curvature Dependent Diffusion Flow on Surface with Thickness 
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Particle diffusion in a two dimensional curved surface embedded in R% is considered. In addition 
to the usual diffusion flow, we find a new flow with an explicit curvature dependence. New diffusion 
equation is obtained in e (thickness of surface) expansion. As an example, the surface of elliptic 
cylinder is considered, and curvature dependent diffusion coefficient is calculated. 
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I. MOTIVATION 

The particle motion on a given curved surface is old but 
interesting problem in wide range of physics. Especially 
the diffusion process of particles on such a manifold is 
still an open problem, and related to various kinds of 
phenomena. 

For example the motion of protein on cell membrane 
has great importance in biophysics. There are several re- 
search papers discussing on this problem. Some of them 
are treating this problem by using usual diffusion equa- 
tion with curved coordinate, and discuss the curvature 
(Gauss curvature) dependence of its solution Other 
of them use the Langevin equation on curved surface and 
calculating the curvature dependence of diffusion coeffi- 
cient @. 

The quantum mechanics of particle motion on such a 
curved manifold is also considered by many authors. This 
problem is usually explained by the Schroedinger equa- 
tion with Laplace-Beltrami operator. However, when we 
treat the curved surface as embedded one in 3 dimen- 
sional Euclidean space, situation is changed and then we 
have a quantum potential term related to the curvature 
additional to the kinetic operator 0, [j] , H . 

Another example is in larger scale physics in which 
our consideration is devoted. Patterns of animal skins 
are well expressed by the reaction diffusion equation @ . 
But the patterns are different for each parts even in one 
individual. For example, Char fish, the side part has 
white spot pattern, but the back part has labyrinth pat- 
tern. (For these two patterns, see for example [7j.) One 
of the reasons might come from the curvature difference 
between side part and back part. If the diffusion is influ- 
enced by the curvature, this difference of patterns might 
be explained. Furthermore, the cross section of fish has 
form of ellipsoid and the surface can be approximated as 
the one of elliptic cylinder. In two dimensional space, 
we have only two kinds of curvature, one is Gauss curva- 
ture and other is mean curvature. Both are constructed 
from second fundamental tensor by taking determinant 
or trace. Gauss curvature can also be constructed only by 



metric tensor and its derivatives, but this is not the case 
for the mean curvature. The elliptic cylinder, in which we 
have much interest, has zero Gauss curvature and non- 
zero mean curvature. Therefore to explain the pattern 
change of Char fish, solution of the diffusion equation 
should depend on mean curvature. This is impossible if 
we start from usual diffusion equation because it depends 
only on metric but not on second fundamental tensor. 
Therefore we need some new diffusion equation, which 
bring not only Gauss curvature but also mean curvature. 
In this article, we discuss how to construct such curvature 
dependent diffusion equation. 



II. COORDINATE AND METRIC 

The simple extension of diffusion equation in Euclidean 
space to Riemannian space can be done by changing 
Laplacian with Cartesian coordinate to the one with 
curved coordinate, i.e. Laplace Beltrami operator. This 
coordinate change is not enough for our purpose, how- 
ever. The way of construction of new diffusion equation 
in this paper is the followings. We re-identify the two 
dimensional diffusion as the limiting process from three 
dimensional diffusion. We place the curved surface E in 
three dimensional Euclidean space i?3, and we put two 
similar copies of E, called E and E' at a small distance 
of e/2. Our particles can only move between these two 
surfaces, and later we take a limit e — * 0. We look for the 
form of diffusion equation in this limit. The coordinates 
we use hereafter is the followings. (See fi.g.1) 

X is the Cartesian coordinate in i?3 . x is the Cartesian 
coordinate which specifies only the points on E. q l is the 
curved coordinate on E. (Small Latin indices k, • • • 
runs from 1 to 2.) q° is the coordinate in R3 normal to 
E. Further by using the normal unit vector ^(g 1 ,^ 2 ) on 
E at point (g 1 ,^ 2 ), we can identify any points between 
two surfaces E' and E by the following thin-layer approx- 
imation [5|. 



X(q°,q\q 2 ) = x{q\q 2 ) + q rt(q\q 2 ), 



(1) 
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where -e/2 < q° < e/2. 

From this relation we can obtain the curvilinear coor- 
dinate system between two surfaces (c R3) by the co- 
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FIG. 1: Embedding and Coordinate 

ordinate = (q° , q 1 , q 2 ), and metric G M „. (Hereafter 
Greek indices /U, u, ■ ■ ■ runs from to 2.) 



_ dX dX 



(2) 



Each part of G M1/ is the following. 



Gij = gij+q"( 



where 



dx dri dx dri 
dq % dqi dqi dq l 



dri dri 
dq l dqi 



dx dx 
9ij = dq 1 'W 



, (3) 



(4) 



is the metric on E. Hereafter indices i,j,k--- are lowered 
or rised by g^ and its inverse g lJ . We also obtain 



Goi = Gm = 0, Goo = 1- 



(5) 



We can proceed the calculation by using the new vari- 
ables. We first define the tangential vector to E by 



Bk = 



dx 
dq k 



(6) 



Note that ri ■ B k = 0. Then we obtain two relations. 
Gauss equation: 



dB 



Weingarten equation: 

dri 



(7) 



(8) 



where 



r« = ^g km (d igmj + d 2 g im - d m g l0 ). 

is called Euler-Schauten tensor, or second fundamen- 
tal tensor defined as 



K-- — -B- 

13 ~ dq* r 



(9) 



The second fundamental tensor t^ij IS the projection of 
dri into the surface. Furthermore, the mean curvature is 
given by 

K = g ij K ij , (10) 
and Ricci scalar curvature R is obtained by 

R/2 = det(g tk Kkj ) = det(«*.) = ^(k 2 - K ijK ij ). (11) 

Then we have the formula for metric of curvilinear co- 
ordinate in a neighborhood of E. 

Gij = gi j + 2q°Kij + (q ) 2 K miK f. (12) 

Under the inversion q° — > — q°, we have 
well as ri — + — ri from ri = 8qX / \ d^X |. Therefore Gij 
is invariant under q° — » — g°. 



Now we have the total metric tensor such as, 

1 



G M „ — 



G 



(13) 



III. EMBEDDING OF DIFFUSION FIELD 

Let us denote 3 dimensional diffusion field as (f>^ , and 
Laplacian as A^ 3 '. Then we have the equation with nor- 
malization condition 



d^ 

dt 



DA< 3 V (3) , 



i = / ^H q °, q \ q 2 )VGd\ 



(14) 
(15) 



where D is the diffusion constant, and G = det(G„ J/ ) = 
det(Gy). Our aim is to construct the effective two di- 
mensional diffusion equation from 3D equation above. 



d^ 

dt 



1 = J^HqW^Jq, 



(16) 
(17) 



where is the two dimensional diffusion field, g = 
det(<7y), and A^ 6 -^ is unknown effective 2D diffusion 
operator which might not be equal to simple 2D Laplace 
Beltrami operator. 

From two normalization conditions, we obtain 

1 = J ^\q°,q\q 2 )VGd\, 

= J[f' 2 dq (^VGT9)]Vgd 2 q, 



r e/2 
-e/2 

4> {2) {q\q 2 )V9 d 2 q . 
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Therefore we obtain the relation, 

e/2 



(2 V,g 2 )=/ t i3) dq°, (18) 

-e/2 



where 



(19) 



We multiply y/G/g to equation (TT4")) and integrate by 
q°, then we obtain 

" -D A^^dg , (20) 



Of 



from (|18[) . we need to expand A< 2 ) up to C(e 2 ). The 
following relations are useful 

Gij = gij + 2q\ tJ + {q°) 2 K imK f, (28) 
G 13 = g ij -2q°K ij + 3( g °) 2 <^ mj + 0(e 3 ), (29) 
G - .9 {1 + 2q°K + (q°) 2 (K 2 + R) + 0(e 3 )}, (30) 

G 1 ' 2 = g^ 2 {l + q°K + l(q°) 2 R + 0(e 3 )}, (31) 

where R = k 2 — KijK 13 is used. 

Then the operator A^ 2 ) can be expanded as follows 



e/2 



where 



A(2) = A( 2 > +q°A + (q ) 2 B + O(e 3 ), (32) 



A^ = y/Gfg A (3) ^?G. (21) 



where, 



Next we analyze new operator A^ 3 ) . From the form of 
Laplace Beltrami operator 

A (3) = G _1/2 — G 1/2 G^ — , 



we have 



= A< 2 >+A« 



(22) 



where 



and 



1 ,,, <9(k 2 - -R) _ j, 9k 
+ TTfl o — - + 2K y — ). (34) 

Then our two dimensional effective diffusion equation 
up to 0(e) is, 

= DA^S {2) 



A< 2 > EE . g -V2|_ G l/2 Gy ^_ (g/G) l/ 2j (23) 



0t 



A (1) = — G 1/2 — G" 1/2 (241 
~ dq° dq° • [ ' 

Then our diffusion equation has form 



dt 



= D A< W« 

J-e/2 



(25) 



f e/2 

+ DA q°4>Wdq° 
J -e/2 
r e/2 

+ DB / (q ) 2 ^dq + 0(e 3 ). (35) 

J -e/2 

To proceed the q° integration, we suppose there is no 
diffusion flow in normal direction in layer , that is, 

^(3) 1/2 dG-^ 

° = ^r = 9 dq° ■ (36) 



Solution is, 



The contribution from A^ 1 ) vanishes because 
f /2 A^4> {3) dq a = g- 1 ' 2 f' 2 A(G) 1 / 2 dq° 

J -e/2 J -e/2 Oq U dq" 

= 9- 1/2 i (G) 1/2 ^l 11%= 0. (26) 

The last equality is the requirement that diffusion flow 
does not pass through the surface: £' and E. 

Now we calculate r.h.s of ([25|) up to 0(e 2 ). Since we 
have 



4> {3) = ^(G/gy/^Hq 1 ^ 2 ), (37) 



V : / 72 (G/g)^ 2 dq°. 

-e/2 



(38) 



^ 3 ) = o( e - x ), 



(27) 



Each integration can be explicitly performed, and we 
obtain 

N = e+^e 3 + 0(e 5 ), (39) 

(40) 
(41) 



<1°> = ^+0{e% 
<(q°) 2 > = T2+ 0{ ^' 
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where we have used the definition 

< f(q°) >^ 1 f' 2 f(q )(G/gy/ 2 dq . (42) 

iV J-e/2 

We obtain the final form of equation up to 0(e 2 ) as 



d^ 2 

at 



= DA^^ +D(Ak + B)<I>W 

= DA^^+Dg- 1 / 2 ^ g 1 / 2 
oq % 

pi i n p 

x {(3« m - 2«««)^ - -g«^ J }^ 2 ), (43) 



where D — j^D- 

We give two comments here. First, 0(e) term disap- 
pears. Since e has the dimension of length, it always 
appears with curvature n. Therefore the 1st order term, 
if it exists, it contains 1st order of curvature k. But this 
curvature depends on unphysical choice of normal unit 
vector n, and so it does not appear. 

Second, additional potential term disappears. In quan- 
tum mechanics, the similar embedding techniques leads 
to the appearance of additional potential term written 
by curvature. But in our classical case we have no such 
terms. Because in diffusion equation, potential term 
breaks probability conservation law, i.e. 

d<j>/dt= (DA + V{x))<p, 



IV. PROPERTIES OF CURVATURE 
DEPENDENT FLOW 

The anomalous flow equals to zero for the flat surface. 
The last term in equation (|45[) shows that curvature gra- 
dient generate the flow without particle density gradient. 
From the signature of this term, this flow goes from the 
smaller Ricci scalar point to the larger Ricci scalar point. 
Ricci scalar can take the negative, zero, and positive val- 
ues. (Ricci scalar R is related to Gauss curvature by 
R/2 = det[/t*-].) Let us work with the coordinate which 
satisfies 



9i. 



Sij, k) = diag[l/ri, l/r 2 ], 



at point P, where r*j is the curvature radius along the 
q l coordinate and it takes positive or negative value for 
convex or concave. (The metric can be diagonalized by 
choosing the two coordinates as to satisfy orthogonality, 
and it can be normalized by using the re-parametrization. 
The second fundamental tensor is diagonalized by rota- 
tion of coordinate.) 



Then at point P, we have R 



and, 



• R > if the surface is convex or concave. 

• R = if the surface is essentially flat. 

• R < if the surface is hyperbolic. 

Therefore the flow goes from hyperbolic or flat points 
to convex or concave points with positive larger Ricci 
scalar value. 



d 3 x cj> = / d 3 xV(x)(j) ^ 0. 



dt 



The normal diffusion flow can be written in general 
coordinate, 



Jn = -Dg l 



(44) 



The anomalous diffusion flow is 
J A = -£>{(3k™<-2k^) 



W 2 ' IA>}.(45) 



The Diffusion equation is written as 



dt 



V,(A + J\), 
■9~ 1/2 #j3 1/2 (^ + JA), (46) 



where V, is the covariant derivative. By using a suitable 
boundary condition, we can prove 



d_ 

It 



J w 



/2 d 2 q = 0. 



Next we consider the first term in (|45|) . We have posi- 
tive or negative value for 

f ij = 3K im Ki n -2KK ij , 

depending on the value of curvature. In our coordinate, 
we can immediately write it in the simple form 



/*i = ( fi ( J- _ JL). 

When the surface is hyperbolic (R < 0), 



li 



1 



(ri) 2 | rir 2 | 



2 >o, r= 1 



(r 2 ) 2 | nr 2 | 



usual diffusion occurs (See fig. 2). 

When the surface is convex or concave (R > 0), 



1 2 = | r 2 | -2 | n 

(ri) 2 | rir 2 | | r\ | 2 | r 2 | 



1-22 



''i 



-2 | 7- 2 

2 kil 



(47) 



>o, 



FIG. 2: Wave packet on hyperbolic surface diffuses in two 
directions. 



FIG. 4: Wave packet on convex surface with one direction 
curvature is over two times higher than another. Packet dif- 
fuses in higher curvature direction and concentrates in smaller 
curvature direction. 



FIG. 3: Wave packet on convex surface concentrates. 




In this case, we have three possibilities. One possibility 
is that both are negative, if 

1/2 <| ^ |< 2. 

n 

Then we have no diffusion but concentration occurs (See 
fig. 3). 

The second possibility is that one of two is positive and 
the other is negative, if 

| - |< 1/2, or | ^ |> 2. 
n r\ 

Then we have diffusion in one direction, but concentra- 
tion in another direction (See fig. 4). 
The third possibility for R > is, 



— = I/ 2 : OT — 

n n 



2. 



This is critical point, where the flow stops for larger 
curvature direction and flow concentrates for smaller 
curvature direction. 

When the Ricci scalar is zero (R = 0), for example 
r 2 = oo, f 22 = and f 11 > follows. The diffusion 
occurs only in q 1 direction but not in another direction. 
(See fig. 5). 

In this way, this anomalous diffusion flow has much 
varieties depending on the curvature. 



ONE EXAMPLE: ELLIPTIC CYLINDER 



FIG. 5: Wave packet on elliptic cylinder. Packet diffuses only 
in curved direction but not in another direction. 



cylinder is the case of R = just as figure 5, but the 
surface has non zero mean curvature. 
Ellipsoid is given by the equation 



(-) 2 + (f) 2 
a b 



1. 



(48) 



Then any points on cylinder are specified by curved 
coordinate 9 and z\ 



x = a cos 9, y = b sin 9. 
Another choice of coordinate instead of 9 is, 
du = yjdx 2 + dy 2 = f(9) d9. 
where f(9) is defined as 



f(9) = Va 2 sin 2 9 + b 2 cos 2 9. 



(49) 



(50) 



(51) 



AC 



Let us consider one simple example where diffusion 
coefficient depends on curvature. The surface of elliptic 



FIG. 6: Elliptic cylinder 
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The length of u is given by 

u(<j>) = [%(9)d6 
Jo 

= b[ Vl-k 2 sin 2 9d0 = bE(k,4>), (52) 
Jo 



T = 0.01 



with k = y/1 - {a/b) 2 , a < b. E(k,<f>) is the Elliptic 
integral of the second kind. (See appendix) 

The total length of u is given by 

U = 4bE(k,ir/2). 

We use the normalized value for u hereafter such that, 

u = u/U, 0<u<l. (53) 

The normal unit vector n is given as 
~ = , x y 

«V( i2 /« 4 ) + (W Py/WfiF) + (y 2 / b4 ) 

(&cos6»,asin6>). (54) 



f(0) 



— = -(-b sin 9, a cos 9) — n. (55) 

90 f f 



— ' OX 

B 9 = — = (-a sin 9, b cos 0). (56) 

O0 

Then we obtain the second fundamental tensor. 

n 9 ™ ab frn\ 

Kee = B e - — = —. (57) 

Then we collect all the necessary quantities as follows 

gee = f 2 , g zz = 1, gez = 0, 

ab ab 
Kee = Y> k zz = k 9z = 0, K =~p- ( 58 ) 

Then we obtain the total diffusion equation expressed 
by the parameters 9 and z. 

where the effective diffusion coefficient depends on mean 
curvature. 

where e = e/(2\/3b). 

Under this equation, we obtain the following particle 
number conservation law. 

j t J dzJ 2 Jd9f{9) 4>W(0,z) =0 




FIG. 7: Snap shot of diffusion process starting from the 
wave packet <j> — sin (iru) as the initial condition. At two 
points (0.25 and 0.75), diffusion occurs quickly and its slope 
is smaller than other. 




0.Z 0.4 0.S 0.3 



T U 



FIG. 8: Mean curvature as a function of u when b/a — 2 



with suitable Neumann boundary condition. 

By using the variable u instead of 9 we have simple 
dimcnsionless equation, 



8t 



d_ 

' du ' 



(60) 



where r = tD/U 2 , r] = z/U, V = e 2 (bn) 
AbE{k, tt/2). 



U 



By using the approximation of elliptic function given in 
appendix, curvature dependent potential V can be writ- 
ten as function of u. The simulation can be done as 
usual diffusion equation. For < u < 1 and < tj < 4 
using periodic boundary condition, since the length of 
f] is larger than one of u, the diffusion in u direction 
occurs fastly and then the diffusion in rj direction fol- 
lows slowly just like <f) ~ a + ^2 k b(k)e~ k T cos(krj). The 
it-directional diffusion can not occur uniformly, because 
at u = 0.25, and 0.75 the diffusive coefficient is higher 
than other points. Therefore the slope of diffusion field 
is small especially at these two points during the diffusion 
process, (fig. 7) 



VI. CONCLUSION 

We have discussed on the diffusion equation on curved 
surface embedded in i? 3 . We obtained the new diffusion 
equation up to 0(e 2 ) which includes anomalous diffusive 
flow additional to the usual one. This anomalous flow 
depends on the second fundamental tensor, and it has 
not only diffusion but also concentration properties de- 
pending on the curvature of its manifold. 

At the point with negative Ricci scalar R < 0, sur- 
face is hyperbolic, and diffusion in both direction occurs. 
(fig-2) 
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When Ricci scalar is positive R > 0, we have three 
possibilities, appearing below is curvature radius in 
each direction (i = 1,2). 

• Concentration in both direction (fig. 3), when 

1/2 <| ^ |< 2. 

7*1 

• Concentration in smaller curvature direction, and 
diffusion in higher curvature direction (fig. 4), when 

|^|<l/2, or |^|>2. 

n n 

• Concentration in smaller curvature direction, and 
no flow in higher curvature direction, when 

|-|=l/2, or |^|=2. 

n n 

When Ricci scalar is zero R — 0, surface is essentially 
flat, but we have finite curvature radius in one direction. 
Then we have diffusion only in this direction (fig. 5). 

In the case of surface of elliptic cylinder, we gave a 
concrete form of equation and we showed the curvature 
dependent diffusion coefficient. 



D u = D(l 



12 



-), D Z = D, 



where k is the mean curvature. In this case curvature 
dependence is simply included into diffusion coefficient. 
However this is not true in general case, where situation 
is much more complicated, and this can be seen from the 
form of anomalous flow. 

The application to the pattern formation by reaction 
diffusion using this obtained equation is not yet finished. 
This will be done in further publication. 



VII. APPENDIX 



We approximate the elliptic integral of the second kind. 



E{k,<j>) = I Vl - k 
Jo 



2 sin 2 9 d9 



(61) 



with 



k = y/l- (a/by 



Under the expansion in powers of k 2 , we obtain the 
power series of Elliptic integral of the second kind. 



E(k, 



n=l 



Since the integration part can be expanded by <f> and 
sin 2ncj), we have 



E(k 7 (j>) — aQ(f> + a n sin2n0. 



with the relation 



a = -E(k,ir/2), 



(63) 



(64) 



a n = (-1)' 



,2E(k,n/2) 



tt/2 



sin 2n(f> E(k, <£)d<i>. (n > 1) (65) 



For the real Char fishes, b/a takes values 1.5 ~ 2.5. 
Then the value of k takes 0.75 - 0.92. When b/a = 2, 
each values of a n is given numerically 

a = 0.771, ai = 0.123, a 2 = -0.00506, a 3 = 0.000558. 

Now we have for u given in l[55 



u/b = E(k, 4>) = aQ(f> + a\ sm2(j) + a 2 sin4(/> + • • • . (66) 
And we rewrite it by using normalized u, 



Ol 0,2 , 

4> = 2-ku sin 2<p sin 40 , (67) 

a a 

where u = u/(AbE(k, 7r/2)). 

The iteration method up to order (ai/a^) 1 gives 



4> = 2ttu - — sin(47ru). (68) 
Then we take the derivative by u in both hand sides. 



1 1 n 2ai (a ~\\ 
- = —(1 cos(4ttm)), 

J bao ao 
where the following relation is used. 



(69) 



$ = J(o) = 



b 2 < 



Then the mean curvature given by 
function of u. 



is obtained as 



bn = 



ab 1 , 2ai , %% q , 

-IT = 757-^3 1 - — cos 47TW 3 , 70 



where (3 — b/a. This function is shown in figure 8. 
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